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Some of the first year of this project was spent in starting it up. A Postdoc, Dr. G^rp 
Stenchikov, was hired and began working on September 1, 1992. A graduate student, Yuhe 
Liu, has been supported on this project for the entire time, and another student, Jianpmg Mao, 
has been working on related studies. 

We have studied the observed effects of volcanic eruptions and published one paper 
(Robock and Mao, 1992; copy attached). We have also completed an experiment in conjunc- 
tion with scientists at the Max Planck Institute in Hamburg, Germany (Graf €t uLy 
et al., 1993; copy attached). Both these papers explain the winter warmning produced by 
large volcanic eruptions. We have submitted another paper to the Journal of Climate (Robwk 
and Liu, 1993; copy attached), diagnosing the climatic effects of GCM runs done at NASA 
GISS. We have completed several GCM runs at Livermore of the NCAR CCMl forced by 
observed SSTs and by observed SSTs and volcanic aerosols and are palyzing the results, and 
conducting more runs. These experiments are done in conjunction with Dr. Karl Taylor there. 
Dr. Stenchikov is conducting experiments to help define the appropriate aerosol response 
using his radiative and radiative-convective models. These experiments are continuing. 

Graf, Hans-F., Ingo Kirchner, Alan Robock and Ingrid Schult, 1992: Pinatubo eruption winter 
climate effects: model versus observations. Report No. 94, Max-Planck-Institut fur 
Meteorologie, Hamburg, Germany, 24 pp. 

Graf, H.-F., I. Kirchner, A. Robock and I. Schult, 1993: Pinatubo eruption winter climate 
effects: model versus observations. Submitted to Climate Dynamics. 

Robock, Alan and Yuhe Liu, 1993: The volcanic signal in GISS 3-D model simulations. Sub- 
mitted to J. Climate. 

Robock, Alan and Jianping Mao, 1992: Winter warming from large volcanic eruptions. Geo- 
phys. Res. Lett., in press. 
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WINTER WARMING FROM LARGE VOLCANIC ERUPTIONS 
Alan Robock and Jianping Mao 
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Abstract. An examination of the Northern Hemi- 
sphere winter surface temperature patterns after the 12 
largest volcanic eruptions from 1883-1992 shows warming 
over Eurasia and North America and cooling over the Mid- 
dle East which are significant at the 95% level. TWs pat- 
tern is found in the first winter after tropical eruptions, in 
the first or second winter after midlatitude eruptions, and 
in the second winter after high latitude eruptions. The 
effects are independent of the hemisphere of the volcant^s. 
An enhanced zonal wind driven by heating of the tropical 
stratosphere by the volcanic aerosols is responsible for the 
regions of warming, while the cooling is caused by block- 
ing of incoming sunlight. 

^ Introduction 

The eruption of Mt. Pinatubo on June 15, 1991, was 
possibly the most sulfur-rich of this century, putting about 
20 Mt of SO 2 into the stratosphere [Bluth et al., 1991]. 
Although conventional wisdom holds that volcanic aerosols 
produce cooling at the surface due to reduction of incom- 
ing solar radiation (and indeed that effect is to be expected 
on a global average on a two- to three-year time scale 
[Robock, 1991]), the winter of 1991-92 was extremely 
warm over North American and Eurasian midlatitudes, and 
very cold in the Middle East. Here we present evidence 
that this pattern of warming and cooling was due to ±e 
volcano. An indirect dynamical mechanism in which 
stronger zonal winds are produced by an enhanced pole- 
to-equator temperature gradient in the lower stratosphere 
caused by aerosol heating is responsible. 

It was first suggested by Groisman [1985] that warm 
winters over central Russia were a consequence of large 
volcanic eruptions. Using surface air temperature data for 
stations in Europe (including the European part of Russia) 
and northeastern North America for averages of 2 or 3 
winters after the volcanic years of 1815, 1822, 1831, 1835, 
1872, 1883, 1902, 1912, and 1963, he showed a significant 
warming over the central European part of Rupia, ^d 
insignificant changes in the other regions, including 


cooling over northeastern North America. In a recent 
update, Groisman [1992] showed the winter patterns after 
the 1982 El Chichon and 1991 Pinatubo eruptions and 
found warming over Russia again, but found wanning over 
northeastern North America. He explained the wa rmin g 
over Russia as due to enhanced zonal winds bringing warm 
maritime air from the north Atlantic over the continent, 
and that these winds were due to a larger pole-to-equator 
temperature gradient caused by polar cooling resulting 
from the eruptions. 

While these studies are intriguing, they raise several 
questions. If the larger pole-to-equator temperature gra- 
dient was caused by polar cooling, why are the tempera- 
tures warm in the polar region over Eurasia, not cool? 
Why is northeastern North America cool after some erup- 
tions and warm after others? What are the volcanic pat- 
terns over the rest of the world? How much of the signal 
found is due to volcanoes and how much to El 
Nino/Southem Oscillation (ENSO) events? 

Lough and Fritts [1987], using tree-ring data and 
selecting 24 volcanic years between 1602 and 1900, found 
winter warming over western North America for the aver- 
age of years 0 to 2 after the eruption years. The pattern 
they found is quite similar to that during ENSO years, and 
they made no attempt to correct for this factor. They also 
found spring and summer cooling in the central United 
States and summer warming over the United States west 
coast. They did not use any volcanoes south of 10®S, 
believing them not important for Northern Hemisphere 
(NH) temperature variations. 

Graf [1992] and Graf et al. [1992] presented results 
from a perpetual January general circulation model (GCM) 
calculation of the effects of stratospheric aerosols on cli- 
mate that showed winter warming over both Eurasia and 
North America, and cooling over the Middle East, Africa, 
and the United States. The forcing used in this GCM 
experiment was a stratospheric aerosol loading in the pat- 
tern of the El Chichon volcano. The mechanism producing 
the warming was an enhanced polar vortex caused by 
warming in the tropical stratosphere of approximately 5®C 
from absorption of solar radiation by the volcanic aerosol. 
The stratospheric wind anomaly propagated down into the 
troposphere producing the patterns found. In the lower 
latitudes, the direct effect of reducing the solar flux 
resulted in cooling. However, in the high latitudes, due to 
the weak solar insolation, the dynamical effect dominated 
producing warming. 



In an attempt to comprehensively examine the winter 
surface temperature response to volcanoes, and to explain 
and reconcile the above results, this study looks at the 
response of NH- winter- average surface temperature follow- 
ing each large volcanic eruption since 1883. We examined 
the global surface temperature patterns and corrected for 
the ENSO signal. We chose to start the study in this year 
because before this time the global surface temperature 
data set has such a sparse distribution that it is too difficult 
to discern a pattern. We first explain how we chose the 
volcanoes, and then describe the surface temperature data 
and how they were analyzed. 

Volcanoes 

Volcanic eruptions which are rich in sulfur gases pro- 
duce long-lasting stratospheric aerosol layers which cause 
global-average stratospheric warming and surface cooling 
for several years. These effects can be regional for a few 
months following the eruption while the aerosol cloud is 
inhomogeneous, but in less than a year the cloud becomes 
fairly uniformly distributed and produces no discernible 
regional forcing, with one exception. In the winter, the 
winter hemisphere has very little insolation, so the stratos- 
pheric warming is confined to the tropics and the summer 
hemisphere. Hence the only spatially-varying forcing is an 
enhanced pole-to-equator temperature gradient in the 
winter hemisphere. The aerosol cloud from tropical erup- 
tions appears in the tropical stratosphere within months, 
but takes on the order of a year to reach the tropics from 
mid- and high latitude eruptions. 

In order to decide which volcanoes to use, we exam- 
ined the dust veil index (DVI) [Lamb, 1970, 1977, 1983; 
Robock, 1991] and the volcanic explosivity index (VEI) 
[Simkin et al., 1981, Newhall and Self, 1982; Bluth et al., 
1991; Self, personal comm.]. Because these indices both 
have well-known problems [Robock, 1991], we examined 
all available ice core, astronomical, and geological records, 
but they were insufficient to define a clearly better index. 
We therefore took the conventional approach used in 
dozens of other studies [Robock, 1991] and chose vol- 
canoes which had DVT (d.v.i./En,aj^) ^ 250 or VEI ^ 5 
(Table 1). 


Temperature Data 


The best global surface temperature data set has been 
produced by the Climatic Research Unit of the University 
of East Anglia [Jones et ai, 1986a, 1986b, 1986c; Jones, 
1988] and updated to include global surface air tempera- 
ture observations over land and sea surface temperatures 
from ship observations [Jones et al, 1991; Houghton et al.., 
1990]. We used this data set of monthly average, gridded 
surface temperatures, updated through July, 1992, which 
was kindly provided by PIul Jones. The data set consists 
of temperature anomalies with respect to the metui for the 
period 1951-1980. 

In Order to remove the influence of long-term varia- 
tions, including the warming observed from the begiiming 
of the record to the end, we used a high-pass Lanczos filter 
[Duchon, 1979] to remove variations with a time-scale 
greater than 10 years at each grid point. We also wanted 
to remove the ENSO signal, and so, by correlating tem- 
peranure anomalies with the standardized Southern Oscilla- 
tion Index [Chen, 1982; Ropelewski and Jones, 1987], also 
kindly provided by Phil Jones, with a lag of 2 months over 
the ocean and 4 months over the land, we subtracted out 
the portion of the signal that was correlated with ENSO. 
Li this way, we can interpret our results as solely a vol- 
canic signal, although this simple procedure to remove 
ENSO effects may still leave in some non-linear effects. 

Northern Hemisphere Winter Volcanic Signal 

Because of our understanding of how the volcanic 
aerosol should affect NH winter temperattires, we exam- 
ined the first winter after the tropical volcanoes (Krakatau, 
Santa Maria [including the effects of Pelee and Soufiiere 
the same year], Agung, Fuego, El Chichbn, and Mt. Pina- 
tubo) and the second winter after the high latitude vol- 
canoes (Ksudach, Katmai, and Bezymianny) to allow the 
aerosol time to reach the tropical stratosphere. For the 
midlatitude volcanoes (Tarawera, Bandai, and Quizapu) we 
also chose the second winter, except that the pattern for 
Tarawera manifested itself in the first winter, so we chose 
that one. Evidently the aerosol from the Tarawera eruption 
reached the tropics during the first year. Since the results 
are the average of 12 volcanoes, this choice does not affect 
the basic conclusions. Note that we chose volcanoes from 
both the NH and Southern Hemisphere (SH), and only 
sorted them according to absolute value of latitude. 

The resulting NH winter (DJF) surface temperature 
patterns for each volcano are shown in Figure 1 along with 



the average for the 12 volcanoes. In all 12 cases, a region 
of wanning is seen over northern Eurasia and over North 
America. Figure 2 shows the average for the 12 volcanoes 
and the regions that are statistically significandy different 
from the average of the winters for all the other years 
between 1866 and 1992 at the 95% level according to a 
Student’s t-test. The wanning over northern Eurasia and 
Nonh America is significant at this level. There is a large 
region of cooling over northern Africa, the Middle East 
and southern Asia, portions of which are significant at the 
95% level. The North American warming is not an 
unremoved ENSO remnant, since there is no ENSO signal 
in this average in the tropical Pacific Ocean. 

Because the 1947 Hekla (64°N) eruption was not con- 
sidered large (DVI = 70, VEI = 4), and the 1980 Ml Sl 
Helens (46°N) eruption (DVI = 500, VEI = 5) had a small 
stratospheric input, we did not include the subsequent 
winters in our average, but they both had a very strong 
pattern that closely matches the others shown in Figures 1 
and 2. Several of the eruptions that we included, Bandai, 
Fuego, and Bezymianny, are of questionable global 
significance because of their small DVI, VEI, or geological 
evidence, so we repeated the analysis without them and got 
essentially the same results, with significant warming over 
northern Eurasia and North America. 

As a fimher test of the significance of the winter pat- 
tern (Figure 2), we subjectively examined every winter 
from 1866 through 1992. Out of the 127 winters, 39 had 
patterns similar to that in Figure 2, including all the 12 
winters in this paper and 3 of the winters immediately fol- 
lowing these. Thus, this pattern can occur without vol- 
canic eruptions (in 21% of the non- volcanic years), but 
large volcanic eruptions cannot occur without producing 
this pattern. 


Discussion and Conclusions 

The NH winter surface temperature pattern found 
after large volcanic eruptions almost exactly matches the 
pattern found in the perpetual January general circulation 
model (GCM) calculation of the effects of stratospheric 
aerosols on climate [Graf, 1992; Graf et al., 1992]. The 
warming over both Eurasia and North America, and the 
cooling over the Middle East are all seen in the GCM out- 
put. The patterns also match the European pattern found 
by Groisman [1985, 1992] and the North American pattern 
found by Lough and Fritts [1987]. Thus we conclude that 


winter warming over Northern Hemisphere mid- and high 
latimdes and subtropical cooling are indeed volcanic 
effects to be expected in the winter following large tropical 
volcanoes, in the first or second winter following midlati- 
tude volcanoes, and in the second winter following high 
latitude volcanic eruptions from either hemisphere. The 
forcing caused by heating of the tropical stratosphere pro- 
duces anomalously strong zonal winds in the mid- and 
high latitudes which advect warmer maritime air over the 
continents, producing the warming. In the tropics and sub- 
tropics, the direct effects of blocking of solar radiation by 
the aerosols produces cooling. 

The effects of the mid- and high latitude volcanoes do 
not depend on the hemisphere of the volcano. For both 
NH and SH volcanoes, after the aerosol reaches the tro- 
pics, the dynamical feedback produces winter warming 
over NH continents. In the SH, the lack of continents in 
the mid- and high latitudes prevents the warming effect 
from happeiung in the SH winter. 

We examined the other seasons after volcaruc erup- 
tions and found no other significant regional responses. 
The remarkable cooling in eastern and central North Amer- 
ica and warming in the west in the spring and summer of 
1992, while agreeing with the results of Lough and Fritts 
[1987], are evident in our data, but not at a significant 
level. 

The winter warming in certain locations and global- 
average cooling that lasts for several years are characteris- 
tic features of the climatic system response to volcanic 
eruptions that must be considered when interpreting the 
climatic record. In order to understand the response to 
ENSO events and greenhouse warming these patterns must 
be taken into account. 
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Fig. 1. Northern Hemisphere winter (DJF) surface temperature anomalies for the indicated years (the 
winter following the tropical volcanic emptions, the first or second winter following the midlatitude 
eruptions, and two winters following the high-latitude ones), and the average of the 12 cases. 

Fig. 2. Average of the 12 winter surface temperature 
anomalies shown in Fig. 1 (same as last map in Fig. 1, but 
with contour interval = 0.5°C), Avith shading indicating the 
regions for which this average is significantly different 
from the average of the winters for all the other years 
between 1866 and 1992 at the 95% level according to a 
Student’s t-test. 


TABLE 1. List of 12 largest volcanoes since 1883. 


Volcano Eruption Year Lat. 

d.v.i./Emsx 

VEI 

Krakatau 

1883 

6°S 

1000 

6 

Tarawera 

1886 

38°S 

800 

5 

Bandai 

1888 

38°N 

500 

4 

Santa Maria 

1902 

15°N 

600 

6 

Ksudach 

1907 

52°N 

500 

5 

Katmai (Novarupta) 

1912 

58°N 

500 

6 

Quizapu (Cerro Azul) 

1932 

36°S 

70 

5 

Bezymianny 

1956 

56°N 

30 

5 

Agung 

1963 

8°S 

800 

4 

Fuego 

1974 

14°N 

250 

4 

El Chichdn 

1982 

17°N 

800 

5 

Ml Pinambo 

1991 

15°N 

1000 

6 



Krakatau 1883-84 


Tar awe r a 1886-87 


Banda i 1889-90 
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